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ABSTRACT 

The transiting planet CoRoT-lb is thought to belong to the pM-class of planets, in which the thermal emission dominates in the optical wavelengths. 
We present a detection of its secondary eclipse in the CoRoT white channel data, whose response function goes from ~400 to ~1000 nm. We used 
two different filtering approaches, and several methods to evaluate the significance of a detection of the secondary eclipse. We detect a secondary 
eclipse centered within 20 min at the expected times for a circular orbit, with a depth of 0.016±0.006%. The center of the eclipse is translated in 
a 1-cr upper limit to the planet's eccentricity of ecosw <0.014. Under the assumption of a zero Bond Albedo and blackbody emission from the 



planet, it corresponds to a Tc 



=2330*S K. We provide the equilibrium temperatures of the planet as a function of the amount of reflected light. 



If the planet is in thermal equilibrium with the incident flux from the star, our results imply an inefficient transport mechanism of the flux from the 
day to the night sides. 

Key words, planetary systems - techniques: photometric 
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1. Introduction 

When a Hot Jupiter transits its host star, unless the eccentric- 
ity is anomalously large, it will also produce secondary eclipses, 
sometimes also called occultations, from which precious infor- 
mation about the atmospheres of these intriguing objects can 
be inferred. During the last four years, we have witnessed the 
first detections of thermal emission from several Hot Jupiters 
at diff erent band passes in the infrared, most of them from 
space (Charb onneau et ai1l2005l: fPeming et al. 2005). These pi- 
oneering studies have revealed some features of their atmo- 
spher es, such as thermal i nversi ons at hight atmospheric alti- 
tudes dKnutson et al.ll2008ll2009l) . 

The different opacity of the planetary atmosphere, especially 
due to the abundances of TiO and VO, can m ake the spec- 
trum emitted by the planet very different, what led lFortnev et alj 
(2008) to propose two classes of exoplanets, pM and pL. Those 
receiving more incident flux from their star, the pM class, would 
exhibit thermal inversions of their stratospheres, and as a con- 
sequence the depths of the secondary eclipses in the optical and 
near infrared bands will be bigger than expected for a pL planet. 
These authors, as well as Lopez-Morales & Seager (2007), argue 
that in the red parts of the optical spectrum, the thermal emission 
from these objects is much more important than the contribu- 
tion of the reflected light, as theoretical mo dels are in agreement 
with a very low albedo A for pM plan ets dSudarskv et alfe OQO; 
iBurrows et all2008HHood et al.ll2008l) . 
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* Based on observations obtained with CoRoT, a space project op- 
erated by the French Space Agency, CNES, with participation of the 
Science Programme of ESA, ESTEC/RSSD, Austria, Belgium, Brazil, 
Germany and Spain. 



The first transitin g planet detected from space, CoRoT- 
lb (|Barge et all 120081) . completes an orbit around its GOV star 
every 1.5 d. Its incident flux puts this planet in the pM cate- 
gory. In this paper, we focus on the detection of its secondary 
eclipse. We describe the data set and the preparation of the light 
curve for the search for the secondary in Section |2l and the dif- 
ferent techniques used to evaluate the significance and depths of 
the secondary in Section [3] Finally, in Section |H we discuss the 
physical implications of our results. 

2. Observations 

CoRoT-lb was observed during the first observing run of 
CoRoT, attaining a total duration of 52.7 d. In the present 
study, the photometry was performed using the latest version of 
the pipeline, which uses the information about the instrument's 
Point Spread Function and the centroids of the stars measured in 
the asteroseismic channel to correct for the effects of the satellite 
jitter in the white light curve. The effect of the Earth eclipses, in 
which a thermal shock is translated in a bigger jitter and thus less 
flux inside the photometric apertures, is greatly improved in this 
new version of the pipeline. 

The light curve was sampled every 512 s for the first 28 days, 
and then changed to 32 s until the end of the run. During the 
512 s part, the dispersion of the normalized data is 0.00079, and 
0.0019 in the 32 s section. The pipeline version of the light curve 
is plo tted in Fig . Q] To evaluate the red noise content in our data 
(e.g. lPont et al.l2006l) . we computed the standard deviation of the 
light curve (filtered with a median filter with a window of 12 h) 
with different bin sizes, ranging from 30 to 500 min. The result is 
plotted in Fig. [2j the 32-s part of the curve shows less noise than 
the 512-s part below bins of ~ 150 min, while it achieves similar 
levels as the 512-s part above that size. This reveals the pres- 
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ence of uncorrected low frequencies at levels of <1Q 4 in units 
of normalized flux, and of uncorrected signals of about 2x1 0~ 4 
with periods between 30 and 150 min in the 512-s part of the 
data. 
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Fig. 1. Light curve of CoRoT-lb prepared for secondary eclipse 
search (the transits have been removed). For displaying pur- 
poses, the data were combined in 10-points bins. The black 
points show the pipeline version of the curve, while the blue 
points show the curve corrected for the orbital residuals, hot pix- 
els and remaining outliers, as described in the text. 



3.1. First analysis: filtering, local polynomial normalization 
and trapezoid fit 

Two jumps in the data, due to hot pixels, were corrected by es- 
timating the median fluxes before and after the jump, and we 
rejected the points belonging to the region between CoRoT date 
(HJD - 2451545) 2637.1 and 2637.75, as they were affected by 
the stabilization of the hot pixel to its normal levels. 
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Fig. 3. Amplitude spectrum of the light curve of CoRoT-Exo-1, 
before (left) and after (right) the correction for the orbital resid- 
uals described in the text. 
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Fig. 2. The standard deviation of the light curve after removing a 
median filtered version with a 12-h window, and after combining 
it in bins of different sizes. Black dots are for the 32-s part of 
the light curve, and blue dots for the 512-s part. The solid line 
follows the theoretical diminution of the dispersion from the 30 
min point in the 32-s part, if the noise were white. 



3. Analysis of the CoRoT light curve 

We present two different analysis of the data in order to search 
and measure the secondary eclipse of CoRoT-lb. They differ 
substantially from the preparation of the light curve to the es- 
timation of the significance of the secondary eclipse detection, 
and we considered both in order to reinforce or invalidate a de- 
tection of the secondary eclipse. 



The signal we want to measure is at the 10 level, and thus 
we inspected the level of noise in the li ght curve by producin g 
its amplitude spectrum using Period04 dLenz & Bregerl 12005). 
In the left panel of Fig. [3] some peaks at the orbital frequency 
and its daily aliases remain at a level of 1.5 x 10~ 4 . This mo- 
tivated us to try to improve the filtering at the orbital period. 
We estimated the signal at each satellite's orbit i from the sig- 
nal in the 30 orbits before and after the orbit i, as described in 
lAlonso et al.1 {2008). Some remaining outliers were located by 
subtracting the low frequencies in the curve using a moving me- 
dian filter with a window of about 1 d, and flagging iteratively the 
points at more than 3 times the standard deviation calculated in a 
robust way. Finally, these points were replaced by interpolations 
to their closest points. In total, these interpolated points repre- 
sent 3.0% of the original data. The amplitude spectrum of the 
corrected curve is plotted in the right panel of Fig. [3] where the 
amplitudes of the peaks at the orbital frequency and its aliases 
and overtones have clearly been reduced. 

To avoid possible effects of unequal number of points at dif- 
ferent orbital phases due to the phases at the beginning and end 
of the observations, we cut the parts of the data before the first 
recorded transit and after the last one. As we are interested in the 
secondary eclipse, we also cut the transits in the curve. The light 
curve prepared for the search for the secondary eclipse is plotted 
in Fig.Q] 

For the search of the secondary ecli pse, we followed the 
method described in I Alonso et al.l (l2009al) . Basically, it consists 
in evaluating the depths of fits to trapezoids with the same over- 
all shape as the planet's transit, at different phases of the or- 
bital period. At each test phase 0,-, the low frequencies are re- 
moved by performing, at each planetary orbit, linear fits to the 
regions [0; - <p a - 4> b ,4>j - <p a ) U(0; + 4>a,4>i + 4>a + 4>b\, where 
<p a and (pt, were chosen in order not to include points inside 
a transit duration, but at the same time allowing enough data 
points to enter in the fit. We obtained a good compromise us- 
ing the values </> a = 0.04 and 0/,=O.12. Once all the planetary 
orbits were normalized this way, we binned the data with a 
size of 0.001 in phase, and fitted t he depth of a trapezoid using 
a Levenberg/Marquardt algorithm dLevenber jdl 1 944t [Marquardt 
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Fig. 4. Depth of a trapezoid with the shape and duration of a 
transit, as a function of the planet's orbital phase. The maximum 
is centered in phase 0.5, corresponding to the secondary eclipse. 
In grey, the result in a light curve where the secondary eclipse 
signal has been diluted (see text for details). 
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Fig. 7. The x 2 f° r the duration of the secondary eclipse, and the 
1,2 and 3cr confidence limits. The vertical line shows the total 
duration of a transit. The detected signal has the same duration 
as the transit at 1-sigma level. 
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Fig. 5. Phase folded curve of CoRoT-lb during the phases of 
secondary eclipse. The data have been binned in 0.01 in phase 
(~22 min). 
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Fig. 6. The x 2 space for different centers and depths of the sec- 
ondary eclipse, and the different <x confidence limits. 



1 1963b . fixing the rest of the trapezoid's parameters to the val- 
ues of the transit. The final depth of the fits as a function of the 
orbital phase is plotted in Fig. 2J where the maximum is well 
centered at the expected orbital phase 0.5. We can evaluate the 
significance of this detection by computing the dispersion of the 
fitted depths in the parts of the phase diagram not affected by 
the inclusion of the secondary eclipse at phase 0.5 in the regions 



where the fits used to normalize the transits were computed, i.e., 
<p e [0, 0.5 - (p a - <pb) U(0-5 + (p a + (f>b, 1] - The significance of the 
detection calculated this way results in 2.1-cr. The phase folded 
light curve around the secondary phase and the best fit trapezoid 
(with duration and shape fixed to that of the transits) are shown 
in Fig. |5] 

We performed the same technique in a light curve where we 
diluted the signal of the secondary eclipse. To do so, we sub- 
tracted from the light curve a version of it, filtered with a mov- 
ing median using a window of ~1 d in duration. We shuffled 
randomly the residuals, and we added back the subtracted fil- 
tered curve. The resulting depth vs. phase diagram using the 
method described above in this curve is plotted as a grey line 
in the Fig. [4] If we take the dispersion of this measured depth 
(3.29xl0~ 5 ) as the precision in the measurement of the sec- 
ondary eclipse, then the significance of the signal at phase 0.5 
is 4.5-cr. 

Additionally, three different methods to evaluate the depth 
and significance of the secondary eclipse signal were tested. In 
the first method, we removed the best fit trapezoid to the data 
(where the only fitted parameter was the depth, the rest was 
fixed to the values of the transit), shifted circularly the residu- 
als, reinserted the signal, and re-evaluated the fitted depth. The 
final depth thus takes into account the effect of red noise in the 
data, and it is of 0.014+0.002%. The second method consisted in 
fitting two gaussians to 1) the distribution of points inside total 
eclipse and 2) a subset of the points outside the eclipse with the 
same number of points as in 1), and compare their fitted centers. 
This fit was performed to 500 subsets of 2), with a randomly cho- 
sen starting point and varying the size of the bins in the distribu- 
tion between 0.005% and 0.02%. The result is of 0.021±0.003%. 
In the third method, we explored the x 2 distribution in a grid 
of centers (from -60 to +60 minutes from the expected center) 
and depths (from 0.01% to 0.04%) of the secondary eclipse. The 
other parameters of the trapezoid were fixed to the values of the 
transit. The minimum ^ 2 and the 1,2 and 3-cr confidence lev- 
els are presented in Fig. [6] and the best fitted depth using this 
method is 0.016+0.006%. We show the x 2 map of the duration 
of the trapezoid and the 1,2, 3-cr confidence levels in Fig. [7] The 
duration of the secondary eclipse is, as expected, compatible at 
a <2-cr level with the duration of the transits. 



3.2. Second analysis: using the IRF filter 

The analysis presented so far shows encouraging evidence for 
the detection of a secondary eclipse at the expected phase and 
duration for a circular or almost circular orbit. However, the de- 
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tection is unavoidably tentative given the extremely shallow na- 
ture of the signal. Each step of the analysis involved a number 
of free parameters, from the hot pixel and satellite orbit resid- 
ual corrections to the individual corrections applied for stellar 
variability local to each putative secondary location. 

In an effort to reinforce or invalidate the detection, we car- 
ried out a separate secondary eclipse search using different pre- 
processing and eclipse detection methods which were designed 
to minimize the number of free parameters. The regions around 
the hot pixel events were simply clipped out, no correction for 
orbital residuals was ap plied, and we used t he ite rative recon- 
struction filter (IRF) of lAlapini & Aigrainl d2009l) . which has 
only two free parameters, to isolate signal at the planet's orbital 
period from other signals including stellar variability. 

The starting point of this analysis was the same light curve, 
as described in Section[2] To circumvent the issue of varying data 
weights associated with different time sampling, the oversam- 
pled section of the light curve was rebinned to 512s sampling. 
Outliers were t hen identified and cl ipped out using a moving me- 
dian filter (see Ai grain et al.ll2009l for details). Finally, we also 
discarded two segments of the light curve, in the CoRoT date 
ranges 2594.25-2594.45 and 2637.10-2637.70), corresponding 
to the two hot pixel events visible in Fig. [TJ This last step was 
necessary as the IRF cannot remove the sharp flux variations as- 
sociated with hot pixels without affecting the transit signal. 

The resulting time-series was then fed into the IRF A fu ll de- 
scription of this filter is given in Alapini & Aigrain (2009), but 
we repeat the basic principles of the method here for complete- 
ness. The IRF treats the light curve {Y(i)} as Y(i) = F(i)A(i) + 
R(i), where {A(i)} represents the signal at the period of the planet, 
which is a multiplicative term applied to the intrinsic stellar flux 
{F(i)}, and {/?(/)} represents observational noise. A first estimate 
of {A(i)} is obtained by folding the light curve at the period of 
the transit and smoothing it using a smoothing length of 0.0006 
in phase, and is divided into the original light curve. This is then 
run through an iterative non-linear filter (Ai grain & Irwini r2004) 
which preserves signal at frequencies longer than 0.5 days to 
give an estimate of the stellar signal {F{i)}, which is assumed 
to be primarily concentrated on relatively long time-scales. This 
signal in turn is removed from the original light curve and the 
process is iterated until the the dispersion of the residuals re- 
mains below 10~ 4 for 3 consecutive iterations, which occurs af- 
ter 3 iterations in this case. 

The free parameters of the IRF are the smoothing lengths 
used when estimating {A(i)} and {F (/)}■ The former represents 
a compromise between reducing the noise and blurring out po- 
tential sharp features associated with the planet, and the latter 
between removing the stellar signal without affecting the plan- 
etary signal. The values chosen here were adopted by trial and 
error and gave the best results when evaluatin g the performace of 
the IR F for transit reconstruction purposes (Alapini & Aigrain, 
120091) . 

We searched the IRF-filtered light curve for secondary 
eclipses using a very simple sliding box algorithm, where the 
model is taken to be 1 outside the eclipse and the in-eclipse level 
is set to the median of the in-eclipse points (the light curve was 
previously normalised by dividing it be the median of the out- 
of-transit points). This is therefore a 2-parameter model, and the 
phase and duration of the eclipse were varied over a 100 x 100 
grid ranging from 0.3 to 0.7 in phase and from 0.02 to 0.2 in 
duration, where the duration is measured in phase units (corre- 
sponding to 0.7 to 7.2 hours). We then evaluate the significance 
of the detection, in arbitrary units, by dividing the eclipse depth 
by o"/ VC^X where sigma is the scatter of the points - outside 
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Fig. 8. Detection significance map. The red thick cross is the sec- 
ondary eclipse duration and phase with the largest detection sig- 
nificance. The blue thin cross is the largest detection significance 
for a secondary eclipse of the same duration as the primary tran- 
sit. 
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Fig. 9. Depth of a box as a function of the planets orbital phase 
(black curve) with the total eclipse duration (0.0697 in phase 
units) providing the highest detection significance. The largest 
significance is found for a secondary eclipse depth of 0.019% at 
an orbital phase of 0.510. The grey dashed lines show the scat- 
ter (+0.006%) of depths measured in a light curve freed from the 
secondary eclipse with shuffled residuals (grey curve). This scat- 
ter is an estimate of the uncertainty on the depth measurements. 



the primary transit - in the IRF-filtered light curve and N is the 
number of point used to calculate the level inside the secondary 
eclipsf]. The resulting confidence map is shown in Fig. [8] where 
the phase and duration giving the best significance are marked 
with a red cross. The best duration corresponds almost exactly 
to the full duration of the transit (which we measure as 0.0688 in 
phase units from a trapezoidal fit), and the best phase is 0.510, 
in good agreement with the results of Fig. [6] The diagonal pat- 
terns in the significance map are due to individual features in the 



1 The significance can be estimated as S=depth/cr x V^V, and in this 
case the dispersion cr, and the total number of points N remain constant. 
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folded, filtered light curve, which influence a wider range of trial 
phases at longer durations. 

Fig.|9]shows the eclipse depth as a function of orbital phase 
at the best duration. The best-fit depth is 0.019%, consistent with 
the results presented above. One should bear in mind when com- 
paring the different results that the eclipse depth was measured 
here relative to the overall median flux rather than to a local es- 
timate of the out-of-eclipse flux. 

To calculate the significance of the detection, we measured 
the depth and associated uncertainty and divided the former by 
the latter. In this case, we calculated the uncertainty as the dis- 
persion (estimated as 1.48xMAD) around the median level of 
the light curve outside the primary transit, and divided this value 
by the square root of the number of points inside the secondary 
eclipse. This results in an uncertainty of 0.005%, i.e. a 3.5 cr 
detection. We also applied the three other methods already pre- 
sented in Section 3.1. All of these give an uncertainty estimate of 
0.006%, which we adopt as the final uncertainty on the eclipse 
depth derived from the IRE 

The fact that we arrived at consistent estimates of the sec- 
ondary eclipse depth and significance using different approaches 
to filter the light curve and evaluate the uncertainty on the depth 
lends confidence to the detection. Our final adopted values for 
the depth is 0.016±0.006%, consistent with all our estimates 
within one sigma. The eclipse center occurs at the expected 
phase for a circular orbit, with an uncertainty of about 20 min. 

4. Discussion 

We have shown that a decrease in the flux of the light curve 
of CoRoT-lb at the phases of secondary eclipse is detected at 
more than 3-cr level. Its duration is, within 1-cr, the same as 
the duration of the transits, and its shallow depth is of only 
0.016±0.006%. We thus interpret this signal as the secondary 
eclipse, detected at the optical part of the spectrum. 

If we were to explain the secondary eclipse detection 
as reflected light, it would imply a geometric albedo of 
A f =0.20±0.08, which is a bigger value than several up- 
per limits provided by other works in differ ent exoplanets 
dCollier Cameron et al.M2002t iRowe et al.1 [2008). Furthermore, 
theoretical models predict a strong optical absorption in the Hot 
Jupiters' atmospheres ( e.g. ISudarskv et al.ll2000t iBurrows et al.l 
2008; iHood et al.ll2008l) . and for the pM-class of planets, ther- 
mal emission dominates by more than an order of magnitude 
the reflected light. Consequently, we suspect that the secondary 
eclipse signature is not produced entirely by reflected light, but 
most probably by a combination of therma l emission and som e 
reflected light, as in the case of CoRoT-2b (lAlonso et al.L l2009b). 

The CoRoT bi-prism allows to recover chromatic informa- 
tion from the signal. In this case, the difference between the sig- 
nificance of the secondary signal in the blue and the red chan- 
nels might clarify if our detected signal is dominated by reflected 
light or by thermal emission in the optical, as we assumed given 
the arguments above. In the first case, the secondary should be 
detected mostly in the blue channel, while in the case of ther- 
mal emission the red channel would carry most of the signal. We 
checked the colors in the CoRoT-1 aperture, but unfortunately 
the residuals of the jitter correction and the noisier individual 
channels did not allow us to conclude on this subject. Q 

2 After submission of the original manuscript of this paper, we be- 
came awar e of an independent analysis of the red CoRoT channel per- 
formed by Snel len et alj |2009), that allowed these authors to detect a 
0.0126±0.0033% secondary eclipse. This result points towards a small 



The equilibrium temperature can be calculated as 
T eq = T*(RJa) l l 2 [f(l-A B )] 1/4 (1) 

which depends on the Bond albedo A B and the re-distribution 
factor / which accounts for the efficiency of the transport of en- 
ergy from the day to the night side of the planet. / can vary 
between 1 /4 for an extremely efficient redistribution (isothermal 
emission at every location of the planet) and higher values for an 
inefficient redistribution, implying big differences in the temper- 
atures at the day/night sides of the planet. 

To translate the measured depth of the secondary eclipse 
into brightness temperatures, we used a theoretical m odel of a 
GOV star from the Pickles library of stellar models dPicklesl 
1998), calibrated in order to produce the same integrated flux 
as a Planck black-body spectrum with a T e ff=5950 K. Under 
the assumption that the planetary emission is well reproduced 
by a black-body spectrum, and as we know the ratio R /; /R+ 
from the transits, we can calculate the temperatures that pro- 
duce secondary eclipse depths compatib le with our resu l t, using 
the CoRoT response function given in Auverg ne et al.1 J2009). 
The brightness temperature calculated this way resulted in 
TcoRoT=2330^[4Q K. If we further assume that the planet is in 
thermal equilibrium and a zero Bond Albedo, this temperature 
favors high values of the re-distribution factor /=0.57 A di- 
agram showing the implied equilibrium temperatures as a func- 
tion of the albedo, where we have assumed for simplicity that 
all the reflected light is detected in the white channel of CoRoT 
(i.e., that A g =Ag, as in Alonso et al. 09b) is shown in Fig.fTOl 
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Fig. 10. The implied equilibrium temperatures from the sec- 
ondary eclipse measurement, as a function of the amount of re- 
flected light. The best fitted depth and the one sigma limits are 
plotted as the solid lines. Different values of the re-distribution 
factors as defined by Eq. Q] are plotted as dotted lines. Zero 
albedo solutions imply high values of the re-distribution factor. 



Non-zero Bond Albedos, where some of the detected flux is 
reflected light from the star, would imply smaller values of /, 



fraction of reflected light in the white curve, as if thermal emission were 
the only source of received flux we would expect a deeper secondary 
eclipse in the red channel than in the white channel. 
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whilst solutions with A>0.30 are not compatible with the mea- 
sured secondary eclipse. All the solutions with equilibrium tem- 
peratures higher than 1900 K are for / values bigger than 1/4. 
Thus, the secondary eclipse detection favors inefficient redistri- 
butions of the incident flux from the day to the night side. 

One may wonder if the k = R p i/R* obtained from the transit 
fit should be different because of the lack of inclusion of the 
emitted light from the planet in the model. Even for an extreme 
case of / = 1/4, the correction to be applied on the k is a factor 
of 3 below the uncertainty given in Barge et al. (2008), and thus 
we do not consider it necessary. 

The measured center of the secondary eclipse, with an un- 
certainty of 20 min, can be used to constrain the e cos u> <0.014 
at a l-o" level. 

Our measured brightness temperature can be used to predict 
eclipse depths in the K s and z bands of 0.25% and 0.014% re- 
spectively. In these bands, ground-based observations have re- 
cently been successful in the detections of s econdary eclipses of 
exopl anets dSing & Lopez-Moralesl 120091: Ide Mooii & Snellen! 
2009), and for the case of CoRoT-lb, the observations might re- 
veal departures from the black-body assumption s uch as the ther- 
mal invers ions observed in several planets (e.g. iKnutson et al.l 
l2008ll2009l) .FI 
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3 After submission of this work. lGillon et alJ (12009) reported a detec- 
tion in the A*,-band of a 0.28^jJ 4 % eclipse, in good agreement with the 
extrapolation of our results. 
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